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1 INTRODUCTION 



The Sunyaev-Zel'dovich effect dSunvaev fc Zel'dovichlll 97CJ1 
is one of the most interesting and promising observational 
tools for cosmology. During the last four decades, many 
works have addressed its usefulness as a probe that can be 
used to determine cosmological parameters (in special when 
combined with other observational diagnostics such as grav- 
itational lensing or X-ray observations), to study the abun- 
dances and evolution of galaxy clusters up to high redshifts, 
to estimate the gas mass fractions inside clusters, to measure 
peculiar velocities of the gas and to map the inner struc- 
ture of galaxy clusters. Some excellent rece nt reviews on the 
phy sics of the SZ effe ct can be found in iRephaelil (I1995T) 
and IfBirkinshawl 1119991). Recent obse rvati onal results have 
been reviewed in iBirkins haw ( 1993) an d ICarlstrom et alJ 
C200Q). Another excellent review that focus in the appli- 
ca tion of SZ effect observations to c osmology can be found 
in ICarlstrom. Holder fc Reesd J2002I) . 

While the theory of the SZ effect physics is very well 
understood, from the observational point of view its study 
is still difficult. This is mainly due to the faintness of the 
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ABSTRACT 

In this work we study the performance of linear multifilters for the estimation of the 
amplitudes of the thermal and kinematic Sunyaev-Zel'dovich effects. We show that 
when both effects are present, estimation of these effects with standard matched multi- 
filters is intrinsically biased. This bias is due to the fact that both signals have basically 
the same spatial profile. We find a new family of multifilters related to the matched 
multifilters that cancel this systematic bias, hence we call them Unbiased Matched 
Multifilters. We test the unbiased matched multifilters and compare them with the 
standard matched multifilters using simulations that reproduce the future Planck mis- 
sion's observations. We find that in the case of the standard matched multifilters the 
systematic bias in the estimation of the kinematic Sunyaev-Zel'dovich effect can be 
very large, even greater than the statistical error bars. Unbiased matched multifilters 
cancel effectively this kind of bias. In concordance with other works in the literature, 
our results indicate that the sensitivity and resolution of Planck will not be enough 
to give reliable estimations of the kinematic Sunyaev-Zel'dovich of individual clusters. 
However, since the estimation with the unbiased matched multifilters is not intrinsi- 
cally biased, it can be possible to use them to statistically study peculiar velocities in 
large scales using large sets of clusters. 

Key words: methods: data analysis - techniques: image processing - galaxies: clus- 
ters - cosmic microwave background 



effect -specially the elusive kinematic SZ effect- and the 
presence of many astrophysical "contaminants" in the fre- 
quencies where the SZ effect is better observed. 

The SZ effect appears as a secondary anisotropy in the 
Cosmic Microwave Background (CMB), that is, as a small 
spectral distortion in the CMB spectrum caused by the 
scattering of CMB photons due to high energy electrons. 
Therefore, it must be disentangled from the CMB fluctu- 
ations themselves, a task that is relatively easy to accom- 
plish in the case of the thermal SZ effect but very difficult in 
the case of the kinematic SZ. Besides, all the contaminants 
(or "foregrounds") that affect CMB observations affect as 
well SZ effect experiments: extragalactic point sources and 
Galactic foregrounds such as dust, synchrotron and free-free 
emissions. As we will show later in this work, even the two 
different types of SZ effects, the thermal and the kinematic, 
can be considered as contaminants one of the other. 

In the last few years several SZ experiments have 
reached the sensitivities and angular resolutions needed to, 
for the first time in history, fully exploit the power of SZ 
observations. Even more experimen ts are under construc- 
tion or planned for the next year s JCaris^oni^^llj2000|; 
ICarlstrom. Holder fc Reesd 120021: lAghanim et alJ bo04ft . 
Ever for these high-sensitivity experiments the observation 



Herranz et al. 



of the SZ effect is not an easy task, and therefore a great 
care must be placed in the analysis and interpretation of the 
data. 

Several different approaches lead to the many different 
SZ data analysis methods proposed in the literature. One 
approach is the component separation, which goal is to ex- 
tract from the data all the different signals that were present 
for their separate study. In this context the SZ effect is just 
one of the products of the separation. Component separation 
techniques typically make use of the different statistical dis- 
tribution and spectral behaviour (i.e. the different frequency 
dependence) of the components. An example of component 
separation method that has been applied to simulations con- 
taining SZ ef f ect is the Maximu m Entropy Method , e.g. 
HoJsgon^^dJ (M^j). As noted bv lDiego et, al.l IJ2002T) and 
Herranz et al.l l|2002a|), component separation methods can 



be very powerful but there is a risk when the main goal is 
to study only one of the sources (in this case the SZ effect): 
an error in the separation of any one of the components is 
easily propagated to the others, and therefore it is neces- 
sary to be extremely careful that all the assumptions made 
(such as the frequency dependence of the components, or 
their statistical independence, etc) are correct. 

Another approach focuses only on the one component 
under study (the SZ in this case), trying to extract it from 
the data or to estimate a few parameters that characterise it. 
For example, in the case of the SZ effect we can try to obtain 
a map of the Compton parameter or just to obtain the po- 
sitions of the clusters, their integrated fluxes, their peculiar 
velocities, etc. We will call this approach detection and/or 
estimation. Detection/estimation techniques consider all the 
other components of the data apart from the one under 
study as interference (noise) and try to minimise its impact. 
Compared with component separation techniques, they tend 
to be more robust in the sense that they do not need 
to model in detail all the different components that form 
part of the "background". Examples of this type of tech- 
niques applied to SZ data are Ba yesian non-parame tric map- 
ping of t he Compton parameter iDiggoet al. 200 j), optimal 
filtering (JHaghnel^fc T egmarkl Il99a : lAghanim et al.lll997t 
iHerranz et al.l2002al) , Markov Chain Monte Carlo (MCMC) 
sampling techniques llHobson fc McLac hlan 200 3) and pa - 
rameter estimation via simulated annealing i Hansenl2004ar) . 

The choice among the above mentioned techniques is a 
delicate issue; all of them have some desirable properties and 
some drawbacks, and depending on the case under study one 
or another (or even a combination of several of them) will be 
more appropriate. The good thing is that they are not mu- 
tually exclusive and in certain cases they can work together 
to obtain better results. In general non-linear methods such 
as MCMC and simulated annealing are more powerful, but 
they are not free of problems. MCMC methods are expen- 
sive and time-demanding from the computational point of 
view. This make them unpractical for larg e blind SZ sur - 
veys. The SASZ simulated annealing code IIHansenl l2004aT) 
is fast but in presence of astrophysical contaminants such 
as extragalactic point sources and Galactic foregrounds it is 
likely to suffer from systematic errors that can be relatively 
large (?). 

Linear "optimal" filters may be not so powerful as non- 
linear detection/estimation techniques, but they have some 
nice properties that may make us consider their use. They 



are simple to understand and implement. They are quite fast 
and not computationally demanding. Moreover, since they 
do not need a detailed modelling of the background statisti- 
cal properties -apart from the fact that it is statistically ho- 
mogeneous and that it can be properly described by second- 
order statistics, namely the power spectrum- they are very 
robust. Moreover, convolutive filters are appropriate to se- 
lect "compact" features (as clusters of galaxies) in the data 
by filtering out the large scales where diffuse components 
manifest and the small scales where pixel noise appear. The 
multifilters (filters that operate simultaneously on all the 
frequency channels of an experiment, taking into account the 
frequency dependence of the SZ effect and the correlations 
among different channels in order to optimise the d etection 
of the SZ effect) presented in lHerranz et al.l i2002al) provide 
a powerful denoising tool (they minimise the background 
contamination) as well as a straightforward amplitude esti- 
mation criterion. In that work two families of multifilters, 
the matched multifilter (MMF) and the scale- adaptive mul- 
tifiltcr (SAMF) were studied and tested on realistic Planck 
simulations. The results showed that both families of mul- 
tifilters lead to unbiased estimators of the central value of 
the Compton parameter of the clusters in the simulations, 
even in presence of realistic contaminants. MMF give higher 
gain factors (the ratio between the variance of the back- 
ground before and after filtering) than the SAMF, keeping 
the same level of spuriou s detections, so they a re preferred 
as our tool of choice. The lHerranz et al] ]2002aT) multifilters 
were designed to detect/estimate the amplitude of the ther- 
mal SZ effect. In the same way they can be constructed to 
detect the more elusive kinematic SZ effect just by changing 
accordingly the frequency dependence of the effect in the 
formulae that lead to the expression of the filter. 

In the previous discussion we have considered the case 
when only the thermal (or kinematic) SZ effect is present 
in the data. But if the two effects appear at the same time, 
will be the multifilter estimator still unbiased? The answer, 
as we will see, is no. Multifilters use information on the 
frequency dependence and the spatial profile of the sources 
at the same time. If two sources have different frequency 
dependence and profile they are likely to be separated very 
well. But if the two sources have the same profile, or the 
same frequency dependence, the problem will be somewhat 
degenerated and the estimator will show systematic effects. 

A case of two sources with the same frequency 
dependence but different spatial profile is evident: the kine- 
matic SZ effect and the CMB spectra are identical (ex- 
cept for second-order relativistic effects). This means that 
some residual contamination from CMB fluctuations will 
remain after filtering, introducing errors in the estimation 
of the kinematic SZ. Fortunately, the CMB fluctuations are 
damped considerably on t he angular sca l es of t ypical galaxy 
clusters (see for example iHu fc White! J1997l0 . Moreover, 
since CMB fluctuations can be symmetrically positive or 
negative, on average the statistical systematic bias due to 
this effect is zero. Considered individually, however, all kine- 
matic SZ measurements by mean of multifilters will be af- 
fected by a non-removable error, and little can be done to 
solve this problem. 

On the other hand, the problem that appears when two 
signals with the same spatial profile and different frequency 
dependence are overlapped can be solved. In our case these 
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two signals are the thermal and the kinematic SZ effects. 
We will show here that their superposition leads inevitably 
to systematic errors in the estimation of the SZ parame- 
ters when we use standard matched multifilters. This errors 
can be very large in the case of the kinematic SZ effect and 
can be very dangerous for automatic blind SZ surveys. This 
paper is devoted to the study and cancellation of this dan- 
gerous type of bias. 

The paper is structured as follows. In section |5] we will 
briefly review the MMF and its properties. In section [3] wc 
will show how the non-removable systematic bias appears in 
the MMF estimator when two signals with the same spatial 
profile and location but different frequency dependence are 
present. Section 2] shows how to design families of multi- 
filters that overcome that problem: the unbiased matched 
multifilters (UMMF). In section Q3 we will test the MMF 
and the UMMF in realistic Planck simulations. Finally, we 
will discuss our results in section |S] 



2 THE STANDARD MATCHED MULTIFILTER 

The stand ard matched m ultifilte r (MMF) was intro- 
du ced bv IHer ranz et all J2002al) and, independently , 
bv lNaselskv et al.1 feOO^T Foilowing iHerranz et al.l J2002af) . 
let us consider a set of N astronomical images correspond- 
ing to observations at different frequencies of a given region 
of the sky. A signal with a known frequency dependence is 
embedded in the data, so that the observations can be de- 
scribed with the following model: 



d u (x) — fuS v (x) + n„{x), v—1,. 



,N. 



(1) 



where the generalised noise n v {x) corresponds to the sum of 
the other emission components in the map. The background 
n v {x) is modelled as a homogeneous and isotropic random 
field with average value {n v (x)) = and cross-power spec- 
trum P^ 1 ^ 2 (q) (q = \q\) defined by 



{n vi (q)nt 2 (q')) = P„ 1V2 {q) S 2 D (q 



(2) 



where n„(q) is the Fourier transform of n v (x) and 8% is the 
2-D Dirac distribution. 

The spatial profile s v (x) is usually written as the pro- 
duct of the source amplitude A and a spatial template t(x). 
In general, the template at each frequency Vi is the result 
of the convolution of the normalised profile of the source 
and the antenna beam at that frequency. Hereinafter, we 
assume for simplicity that the source has spherical symme- 
try, x = \x\, but the methods can be generalised for non- 
symmetric profiles. For simplicity as well we will assume 
that the source is centred at the position x — 0. 

Equation Q can be expressed in a more compact way 
as a vector equation, 



d(x) = AF(x) + n(x) 



(3) 



where the elements of vector F are F v — f v T v - Hereinafter, 
we will employ the usual boldface notation for the N- 
component vectors in frequency whereas the "arrow" no- 
tation x, q will refer to 2-dimensional vectors in real space 
or Fourier space. 

The goal is to obtain a set of N linear filters that, once 
applied to the data, allow us to determine the amplitude A of 
the signal with a minimum interference of the contaminant 



noise. The filters should take advantage on the knowledge of 
the frequency dependence of the signal f , the characteristic 
profile of the source r and the cross-power spectrum of the 
noise, that can be directly estimated from the data. The fil- 
ters are applied to the N frequency channels in order to give 
a quantity -called filtered map- at each position that will be 



used as an estimator of A. Let T(x) — {v vi (x) 



,(x)} 



the vector of N filters we are looking for, then the filtered 
map at the position b is given by 



*-E/-*«*fl.-H> 



(4) 



that can be expressed in Fourier space in a very compact 
way 



(&) = / dqe-^d\q)T{q) 



(5) 



where the superscript t denotes vector transposition, and 
the operation inside the integral is the usual scalar product 
of vectors. The variance of the filtered map can be expressed 

as 



(w 2 (b))-(w(b)) 2 



dq T'PT. 



(6) 



The explicit dependence in q has been removed from the 
previous equation for the sake of simplicity of notation. Un- 
less it is necessary to include it, we will do the same in the 
following. 

The MMF is obtained when two constraints are imposed 
on the filters T: 

(i) That w(0) is an unbiased estimator of the amplitude 
of the source, that is, w(O) = A. 

(ii) That io(0) is an efficient estimator of the amplitude of 
the source, that is, the variance of the w values is minimum. 



As was shown in IHerranz et alJ (|2002aJ), the filters that 
satisfy these conditions are given by: 



Tmmf = a P 1 F, a 



dqF*p- 1 F, 



(7) 



where P : is the inverse matrix of the cross-spectrum P. 
Using eq. © we have that for the MMF 



<T W = I dq T MMF PT mmf = ol 



(8) 



3 SIGNALS WITH THE SAME SPATIAL 

PROFILE AND DIFFERENT FREQUENCY 
DEPENDENCE 

The MMF is designed to produce an unbiased and efficient 
linear estimator of the amplitude of a signal with known 
frequency dependence and spatial profile in multiwavelenght 
observations. But what happens when there are present 
two signals with the same spatial profile and location but 
different frequency dependence? Such is the case, for exam- 
ple, of the thermal and kinetic Sunyaev-Zel'dovich effects. 
Both effects are superimposed with identical spatial distri- 
bution (if we do not take into account intracluster gas mo- 
tions that can distort the shape of the kinematic effect) but 
quite different frequency dependence. 
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3.1 Thermal and kinematic Sunyaev-Zel'dovich 
effects 

The thermal SZ effect is due to the inverse Compton scatter- 
ing of CMB photons by free electrons in the hot intra-cluster 
gas. It produces a shift in the CMB spectrum to higher fre- 
quencies in the direction where the cluster is observed. This 
effect is characterised by the comptonization parameter y c , 



He 



m e c 2 



T e (l)n e (l)dl, 



(9) 



where ks is the Boltzmann constant, ot the Thomson sec- 
tion, m e the electron mass, c the velocity of light, T e and 
n e are the temperature and density of the electrons in the 
gas and I is the distance along the line of sight. Ignoring 
relativistic effects, the frequency dependence of the thermal 
SZ effect is given by 

( AT \ , N „,_, . hv 



(y) = y c /(/>), 



(10) 



J SZt Kb±CMB 

where h is the Planck constant and Tqmb the temperature 
of the CMB, and where 






(11) 



The kinematic SZ effect, on the other hand, is due to 
the Doppler effect that arises because of the radial peculiar 
velocity of the cluster v r along the line of sight. The intensity 
of the kinematic effect, not taking into account second-order 
relativistic corrections and considering T e to be constant 
along the cluster, is 



(") 



-//< 



v r m e c 



(12) 



V T ) SZk ksTe 

and the frequency dependence is (again ignoring relativistic 
corrections) constant in AT/T units. The intensity of the 
kinematic effect is in general significantly smaller than the 
intensity of the thermal effect. 

At high frequencies, the expressions given above are not 
accurate and relativistic corrections should be taken into 
account when dealing with real data. In this work we restrict 
the discussion to the non-relativistic case for simplicity. 

3.2 Bias in the MMF estimator 

Combining the two SZ effects our data model is 

d(a!) = (^P) = y c TO) - Vr(x)] + n(f), (13) 

\ J- / SZt+k 

where V — (w r m e c)/(/csT e ) and the vector F is constructed 
using the profile r and the thermal SZ frequency dependence 
in eqs. (tTCll and {[HI . 

Imagine a MMF designed for the detection of the ther- 
mal SZ effect is applied to the data in 1131 . The filter will be 
the one given in eq. Q using the vector F that corresponds 
to eq. 1131 . It is easy to calculate the average value of the 
filtered map: 



(Wt(8)>= fdqy^F'-Vr^T 



MMFt — 



y c a" 1 / dq F 4 P F - ycVoT 1 / dq r'P _1 F = 



y c ~ycV^, (14) 

a 



where the constant j3 is defined as 



/3 = I dqrP'^F. 



(15) 



According with eq. 11141 , due to the presence of the kine- 
matic SZ effect the estimator Wt(0) is not longer unbiased, 
and the estimation of y c will have a systematic error pro- 
portional to V/3/a. 

The same will occur when a MMF designed for the de- 
tection of the kinematic SZ effect. In that case, the frequency 
dependence to be used is a constant whose value can be in- 
cluded in V and therefore f v = 1. The shape of the MMF is 
then 



TjfMFt = 7 P l T, 1=1 dq^P 1 t, 
The average value of the filtered map is then 
<» fe (0)> = I dqy c \¥ l - Vt x ] T MMFk 



(16) 



J/c 7 



dqF^'W-ycV-r 1 / dgr'P'V 



Vc- 



1 



VcV. 



(17) 



Note that since y c and V appear as a product in eq. 1131 . 
(iOfc(0j) is not a direct estimator of V but of the quantity 
y c V. Dividing by y c and changing the signs we obtain 



K(0)) 

Vc 



= v- 



i 

T 



(18) 



and therefore we find that our estimator is again biased. 

As we will see later in section^] in the case of microwave 
observations of the SZ effect, the bias in the MMF estimator 
of the thermal SZ effect is normally negligible, but in the case 
of the kinematic SZ effect it is expected to be very large, of 
the order of the effect itself or even bigger. Therefore, this 
bias must be taken into account and corrected. 



3.3 Cancelling the bias 

There are several ways in which the bias described above can 
be taken into account. The first option is simply to separate 
the thermal and kinematic effects before estimating their 
amplitudes. This goal can be attained using a convenient 
component separation technique, for exa mple the Maxi- 
mum Entropy Methodj^Iobsoiie^alK 1 999|) , non-parametric 
Bayesian de tection dDieeo et alJl2 002Tl o r MCMC sampling 
techniques JHobson fc McLachlanl I2003T) . The problem is 
that in practise no component separation technique works 
perfectly, and some uncontrolled residuals always remain 
that will still bias the estimation of the amplitude of the 
effects. 

A second option consists in using eqs. 11141 and 1181 to 
subtract the bias factor in each case. This is a perfectly licit 
option in principle, but can be a little difficult to carry out 
in a real case. The argument is as follows: to determine the 
bias in y c we need, according to eq. 1141 . to know the value 
of y c V, whose estimation is biased. To remove the bias in 
y c V we need, according to eq. 1181 . to know the value of 
y c , which is biased as well. Moreover, on an individual basis 
(cluster by cluster) the quantities V and y c will be estimated 
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with a non-zero error. Therefore, trying to cancel the bias in 
this way leads us to a vicious circle. In section[K]we will see 
that for the case of the future Planck satellite observations 
the bias in y c will be small, and the vicious circle could be 
circumvented just by considering the estimated y c as the 
correct one. This would however introduce further errors in 
any case. Besides, in other experiments the bias in y c can be 
not so small as in the Planck case is. Therefore, we do not 
recommend this option to cancel the bias of the MMF. In 
the next section we will describe two novel lineal multifilters 
that automatically cancel the bias in the estimation of the 
thermal and kinematic SZ amplitudes. 



4 UNBIASED MATCHED MULTIFILTER 

Instead of following any of the approaches described in sec- 
tion 13.31 let us see if it is possible to design linear filters 
that are nearly as efficient as the MMF and do not show the 
systematic bias described above. To do it we will make use 
of the knowledge on the two distinct frequency dependences 
of the thermal and kinematic SZ effects. The resulting fil- 
ters will be called unbiased matched multifilters (UMMF). 
Specifically, we search for a set of filters specifically designed 
to give the thermal amplitude y c without bias and a different 
set of filters specifically designed to give the kinematic am- 
plitude y c V. 



4.1 Unbiased Matched Filter for the thermal SZ 
effect 

We construct the filters *$!={4>v 1 , ■ ■ ■ )VVjv} so that the fil- 
tered map w^ satisfies the canonical matched filter condi- 
tions plus another one aimed at cancelling the bias due to 
the presence of the kinematic effect. 

(i) The filters, when applied to the thermal part of eq. 
I|13|l . give, at the position of the source, the value of the 
comptonization parameter, that is, 



(oF** = 1, 



(19) 



(ii) The filters, when applied to the kinematic part of eq. 
I13[l . give, at the position of the source, no contribution to 
the filtered map, that is, 



V* = 0, 



(20) 



(iii) The variance of the filtered map, a Wij , , is minimum 
(efficient estimator). 

The two first two conditions ensure that the filtered map 
at the position of the source is an unbiased estimator of the 
thermal SZ effect. The third one is the condition for effi- 
ciency of the estimator. The filters that satisfy the previous 
three conditions are 

-J/^lp-^F-/^), A = q 7 -/3 2 , (21) 

where the constants a, f3 and 7 are defined in equations Q, 
11511 and 1161 . It is straightforward to see that (w*(0)) = y c 
and that the variance of the filtered map is 



A' 



(22) 



4.2 Unbiased Matched Filter for the kinematic SZ 
effect 

Let us now find the filters <I>={<^„ 1 , . . . , <j>v N } so that the 
filtered map w<s> satisfies the canonical matched filter condi- 
tions plus another one aimed at cancelling the bias due to 
the presence of the thermal effect. 

(i) The filters, when applied to the thermal part of eq. 
1131 . give, at the position of the source, no contribution to 
the filtered map, that is, 



/ dq F** 



0. 



(23) 



(ii) The filters, when applied to the kinematic part of eq. 
11311 . give, at the position of the source, the value of the 
product y c V, that is, 



dgV * = 1, 



(24) 



(iii) The variance of the filtered map, a WiS ,, is minimum 
(efficient estimator). 

The two first two conditions ensure that the filtered 
map at the position of the source is an unbiased estimator 
of the kinematic SZ effect. The third one is the condition 
for efficiency of the estimator. The filters that satisfy the 
previous three conditions are 



-P- 1 (-/3F + ar) 



(25) 



where the constants a, j3, 7 and A are the same as in eq. 
12H . It is straightforward to see that (w$(0)} = y c V and 
that the variance of the filtered map is 



0",„« = -r 



(26) 



5 APPLICATION TO SIMULATED PLANCK 
DATA 

As a test of the previous ideas, we now apply both the stan- 
dard and the unbiased matched multifilters to simulated 
Planck observations in order to estimate the thermal and 
kinematic SZ of test clusters placed on the simulations. Our 
goal is just to show by an example how the biases described 
in section El appear and how the unbiased matched multi- 
filters are able to cancel it. Therefore, to keep the example 
simple and clear we will restrict ourselves to ideal condi- 
tions in which the spatial profile of the clusters is perfectly 
known from the beginning. A full study of the performance 
of the filters in the Planck case, including uncertainties in 
the cluster profile, asymmetric profiles and realistic cluster 
distributions is out of the scope of this work and will be 
addressed in the future. 

We will focus specially on filters designed to extract the 
kinematic SZ effect. Note that due to the extreme faintness 
of the kinematic SZ effect the task of detecting it is go- 
ing to be very difficult. In fact, Planck sensitivities, angular 
resolution and noise levels do not make it the best experi- 
ment to study the kinematic SZ effect. In spite of this, we 
choose the Planck mission as our example scenario for two 
reasons. The first one is that Planck's instrumental specifi- 
cations, noise levels and performance have been thoroughly 
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Table 1. Technical characteristics of the 9 simulated Planck chan- 
nels. Column two lists the FWHM assuming a Gaussian beam. 
Column three shows the pixel size in arcmin adopted in our simu- 
lations. In column four the instrumental noise variance per pixel 
is given in AT/T units (thermodynamic temperature). 



Frequency 


FWHM 


Pixel size 


^noise 


(GHz) 


(arcmin) 


(arcmin) 


(AT/T) 


30 


33.0 


6.0 


l.lxlO -5 


11 


24.0 


6.0 


1.1x10-° 


TO 


13.0 


3.0 


2.2x10-° 


100 


9.2 


3.0 


6.1xl0" 6 


143 


7.1 


1.5 


1.0x10"° 


217 


5.5 


1.5 


1.6x10"° 


353 


5.0 


1.5 


4.9xl0" 5 


515 


5.0 


1.5 


4.9xl0" 4 


857 


5.0 


1.5 


2.2xl0" 2 



studied in the literature, and realistic simulations of Planck 
observations are relatively easy to produce and analyse. The 
second reason is that at Planck's angular resolution most of 
the clusters will appear as point sources in the sky. There- 
fore, their observed profile will be practically equal to the 
well-known beam profile of the Planck detectors. Then, the 
assumption that the spatial profile of the clusters is known 
is not so unreasonable. 



5.1 Simulations 

Even though for this work we are going to use very sim- 
plistic toy clusters and non-relativistic SZ effects, we intend 
to simulate the other astrophysical and instrumental sig- 
nals that constitute the generalised noise n(x) in eq. J3J 
in the most realistic way possible. Our simulations include 
the latest available information about the physical compo- 
nents of the emission (CMB, Galactic foregrounds and ex- 
tragalactic point sources) and the technical specifications 
of the different Planck channels. The simulations were per- 
formed in patches of the sky of 12.8° x 12.8°. TableQshows 
the assumed observational characteristics of the simulated 
maps. 

The CYs for the CM B simulation were generat ed us- 
ing the CMBFAST code llSeliak fc Zaldarriagalll996l) for a 
spatially-flat ACDM Universe with Q m — 0.3 and Oa = 0.7 
(Gaussian realisation). The simulated Galactic emission in- 
cludes four components: thermal dust, spinning dust, free- 
free and synchrotron. The details of the Gala ctic foreground 
simulat ions are identical to those described in lHerranz et alJ 
J2002ar) (see references therein). The ext ragalactic point 
source simulations come from the model of iToffolatti et alJ 
J1998T) for a Gaussian realisation with the same cosmological 
parameters used for the CMB simulation. 

Synthetic clusters were simulated using the spatial 
profile 



+ r v \y/ri + x 2 y/rl + x 1 



(27) 



The previous profile is a modification of the basic multi- 
quadric profile in which r c takes the role of the core radius 
of the cluster and r v is a limiting cut scale that can be as- 
sociated with the virial radius of the clusters. For x « r v 



the profile in eq. 1271 behaves like the classical /3 model 
(with j3 — 1/2) whereas for x >> r v the profile quickly 
drops to zero. The profile in eq. 1271 is continuous, gives a 
good approximation for the typical cluster profile and is well- 
behaved in Fourier space. For this work we used r v — 10r c . 



5.2 Preliminary analysis of the simulations 

Given a set of maps corresponding to the simulations de- 
scribed above for the nine Planck frequencies, it is easy to 
estimate the values of the cross-power matrix P and there- 
fore to calculate the values of integrals a, ft and 7 for this 
case. Doing so, we obtain the values 



a 


= 2.3721 x 10 lu 


P ~- 


= -1.4482 x 10 ! 


7 = 


= 1.3694 x 10 9 


A -- 


= 3.0386 x 10 19 



(28) 

According to these numbers and eqs. 1221 and 1261 we should 
expect that after using the filters ^ and "J> the filtered maps 
should have variances a* ~ 6.7 x 10" 6 and ag> ~ 2.8 x 10" 5 , 
in AT/T units. 

Let us remember that for the case of the classical 
matched nmltifilters the expected variances of the filtered 
maps should be (JMMFt = (l/°0 — 6.5 x 10" 6 and 
OMMFk - (1/7) 1/2 - 2.7 x 10" 5 , respectively. Then the 
addition of the bias cancellation constraint in the formula- 
tion of the UMMF leads to a small increment of the variance 
of the filtered maps, but a very small one indeed. 

On the other hand, from equations 1141 and 1181 it is 
possible to calculate the expected biases that will arise with 
the standard MMF. Let us denote the estimation of a quan- 
tity £ with the standard MMF as £. We have for the case of 
the Planck simulations that 



v_ 

lie 

v 

V 



77 = 1 + 



1 + 0.06V, 
1.05 



V 



(29) 



Depending on the value of V, some of the previous quan- 
tities can be considerable. Typical values are V ^ 0.1, and 
therefore we find that for the case of the thermal effect the 
bias will be negligible (y c /yc — 1) whereas for the kinematic 
effect the bias is huge (V/V — 10) and can not be discarded. 
Standard MMF for the estimation of the kinematic SZ effect 
is strongly biased. 

The expected errors in the estimation of the amplitudes 
of the thermal and kinematic effects by means of the multi- 
filters is directly related to the variances of the filtered map 
calculated above. Note that in the case of Planck the vari- 
ances in the maps filtered with "kinematic-type" multifilters 
(both $ filters and MMFk) are higher than the variances for 
the "thermal-type" multifilters. That is bad news. It implies 
that the detection and estimation of the kinematic effect by 
means of multifilters is difficult not only due to the faintness 
of the effect, but also to the relatively high intrinsic variance 
of the filtered map. 
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MMFt 
=1.5' ;y =0.0001 ; V— 0.1 



r =1.5'; y -0.0001 ;V— 0.1 




Figure 1. Normalised histogram of the values of the estimated 
y c parameter using standard MMF for the thermal SZ effect. The 
true value of y c is shown with a vertical red line. 




Figure 2. Normalised histogram of the values of the estimated 
y c parameter using ^ UMMF for the thermal SZ effect. The true 
value of y c is shown with a vertical red line. 



5.3 Test 

In order to test whether the numbers obtained above are 
correct or not we ran a set of simulations, including syn- 
thetic clusters with the profile 1271 and the frequency law 
11311 . We simulated 5000 of these clusters. In order to avoid 
cluster overlapping only five cluster per simulated sky map 
were included. For this test, we simulated very bright and big 
clusters, with fixed parameters r c = 1.5 arcmin, y c = 10~ 4 
and V = —0.1. For a temperature of the electrons T e ~ 5 
keV that value of V corresponds to a velocity along the line 
of sight of v r ~ 300 kms^ 1 . The simulations were filtered 
with standard matched multifilters designed for the detec- 
tion/estimation of both thermal and kinematic SZ effects as 
well as with the <& and $ unbiased matched multifilters. 

We will focus here in the performance of the filters re- 
garding the estimation of the SZ effects. We will assume 
that the presence of the clusters is already stablished and 
their locations well known. We would like to remark that in 
a more general situation, where detection must be achieved 
before estimating the parameters, other sources of system- 
atic bias different from the one discussed in this work may 
arise due to the detection criterion itself. For example, if de- 
tection is done by looking for local maxima in the images, 
clusters located in areas where the background is positive 
will be favoured for detection, and then in the background 
will not have zero mean in the areas where clusters are de- 
tected, leading to a new source of bi as. For a more complete 
disc ussion on this, see for exam ple lllerranz et al.l l|2002bh 
and lLopez-Caniego et alj (I2004T) . These bias depend on the 
choice of the detection device and their study is out of the 
scope of this work. 



5.3.1 Performance of MMFt 

Figure shows the results of the estimation of y c for the 
5000 simulated clusters using the standard MMF for the 
thermal SZ effect. The average value of the estimated y c is 
(y c ) — 9.77 x 10 -5 and the dispersion around this value is 



= 5.4 x 10" 



The average standard deviation of the 



filtered maps is, for comparison, a u 



6.2 x 10"°. The 



value for this quantity that was predicted in section IB~2l was 
6.5 x 10~ 6 . Eq. J5U predicts y c /y c = 0.994, and the value 
that is found in the simulations is y c /yc = 0.977. 



5.3.2 Performance of 9 multifilters 

Figure shows the results of the estimation of y c for the 
5000 simulated clusters using the unbiased 9 multifilters for 
the thermal SZ effect. The average value of the estimated y c 
is (yc) = 9.83 x 10~ 5 and the dispersion around this value 
is a Vc = 5.7 x 10~ 6 . The average standard deviation of the 
filtered maps is, for comparison, a Wi , — 6.5 x 10 -6 . The value 
for this quantity that was predicted in sect ion \E7Z\ was 6.7 x 
10 -6 . A small bias is present as well in this case, that can 
be attributed to the finite size of the sample. Hence, for this 
case both the 9 and the standard MMFt multifilters perform 
similarly well, down to the precision of the simulations. 



5.3.3 Performance of MMFk 

Figure [3] shows the results of the estimation of V for the 
5000 simulated clusters using the standard MMF for the 
kinematic SZ effect. The average value of the estimated V 
is (V) = 1.05 and the dispersion around this value is o> = 
0.24. The average standard deviation of the filtered maps 
is, for comparison, o- WMMFk /y c = 0.236. The value for this 
quantity that was predicted in section l5~^l was 0.27. Eq. 1291 
predicts V/V — —9.5, and the value that is found in the 
simulations is V/V = —10.5. 

The previous results have been obtained without pre- 
vious knowledge of y c , that is, the estimation of V for each 
simulated cluster is obtained by dividing the filtered map at 
the position of the cluster by the value y c estimated with the 
MMFt. Therefore, the errors in the estimation of y c do prop- 
agate and affect the estimation of V . If we assume that y c 
is perfectly known we can use is nominal value (y c = 0.0001 
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MMFk 
=1.5' ;y =0.0001 ; V— 0.1 




Figure 3. Normalised histogram of the values of the estimated 
V parameter using standard MMF for the kinematic SZ effect. 
The true value of V is shown with a vertical red line. 
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PHI filters 
r c =1.5' ;y c = 0.0001 ;V=-0.1 



Figure 5. Estimation of the V parameter using <E> multifilters 
for different values of the comptonization y c . The true value is 
indicated with an horizontal dotted line. 




Figure 4. Normalised histogram of the values of the estimated 
V parameter using <J> multifilters for the kinematic SZ effect. The 
true value of V is shown with a vertical red line. 



in this case) and then we obtain (V) = 1.03 and cry- — 0.23. 
We conclude that the uncertainty on y c is relatively unim- 
portant in this case. 



5.3.4 Performance 0/$ multifilters 

Figure0]shows the results of the estimation of V for the 5000 
simulated clusters using the $ multifilters for the kinematic 
SZ effect. The average value of the estimated V is (V) = 
—0.08 and the dispersion around this value is ay = 0.26. 
The average standard deviation of the filtered maps is, for 
comparison, <J w ^/y c = 0.249. The value for this quantity 
that was predicted in section 15.21 was 0.28. The bias has 
been very much reduced, and the small deviation from zero 



can in our opinion be attributed to the limited size of the 
sample and the large variance of the filtered maps. 

If we assume that y c is perfectly known we can use 
is nominal value (t/ c = 0.0001 in this case) and then we 
obtain (V) = —0.08 and ov = 0.25. We conclude that the 
uncertainty on y c is relatively unimportant in this case as 
well. 



5.4 Further tests on the $ multifilters 

The toy clusters used above have large values of the comp- 
tonization y c and the peculiar velocity. Only a few clusters 
in the sky are expected to have such large values of these 
parameters at the same time. Though they serve well for 
our academic test, it is necessary to determine if the filters 
are still unbiased in more realistic cases where the cluster 
parameters take lower values. We performed two sets of si- 
mulations varying the value of y c and V . The same as before, 
for each value of y c and V we simulated 5000 clusters, only 
five of them per map. The core radius was fixed to r c = 1 
pixel. 

Figure|^|shows how the V parameter is estimated when 
varying the comptonization y c . The true value of the velocity 
parameter was fixed to V — 0.1. The estimation is unbiased 
even for low values of y c . The error bars grow as y c decreases, 
which is not unexpected since av = a WiS> /y c . 

Figure|S|shows how the V parameter is estimated when 
fixing y c = 10~ 4 and changing the velocity of the cluster. 
The velocity can be positive or negative. The estimation 
is unbiased for all the considered cases, even for the case 
V = 0. 
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Figure 6. Estimation of the V parameter using <E> multifilters for 
different values of the true value of V. The true value is indicated 
with a dotted line. 



6 DISCUSSION 

In this work we have studied the performance of multifilters 
used to enhance the Sunyaev-Zel'dovich effect signal in CMB 
maps such as the future Planck satellite will obtain. The aim 
of multifilters is to boost the cluster signal with respect to 
the background (CMB plus other astrophysical sources plus 
noise) so that they can be more easily detected and studied. 
Multifilters can be used alone or as a previous step before 
more sophisticated data analysis tools are applied to extract 
the maximum possible amount of information from the data. 
A problem arises when two different signals with 
different frequency dependence but the same spatial profile 
appear in a given position. Such is the case of the ther- 
mal and kinematic SZ effects. In that case, we have shown 
that standard matched multifilters are intrinsically biased. 
This bias can be very strong in the case of the weak kine- 
matic SZ effect. We have designed a family of unbiased 
matched multifilters that cancellate this bias. The price paid 
for this bias cancellation is a small decrease on the gain fac- 
tor with respect to the standard matched multifilters. Unbi- 
ased matched multifilters use the a priori knowledge on both 
thermal and kinematic SZ effect's frequency dependences si- 
multaneously in order to optimise the accuracy of the esti- 
mation. We have tested the performance of both standard 
and unbiased matched multifilters using realistic simulations 
of the nine Planck's frequencies. The simulations contain 
CMB, the main Galactic foregrounds, extragalactic point 
sources and instrumental noise with the levels expected for 
the different Planck channels. With these simulations we 
have shown that standard matched multifilters can estimate 
the kinematic SZ effect with systematic biases much larger 
than the true value of the effect. This problem is automati- 
cally solved using unbiased matched multifilters. The results 



of the numerical tests agree very well with the theoretical 
expectations. 

Our tests show that the statistical error bars in the de- 
termination of the cluster parameters y c and V (or, equiv- 
alently, the cluster velocity v r ) are directly related to the 
standard deviation of the filtered map; in the case of the 
thermal SZ effect this last quantity does not depend on the 
intensity of the SZ effect but on the properties of the back- 
ground (i.e. the cross-power spectrum matrix) and the shape 
of the filter (given by the source and beam profiles) . For the 
Planck sky patches here considered and using ^ multifil- 
ters, the typical error bars in the determination of y c are 
a vc ~ 6 x 10~ 6 , independently of the value of the comp- 
tonization of the clusters, which means that for "bright" 
clusters (y c ~ 10 -4 ) individual comptonizations can be de- 
termined with errors ~ 20% using the Planck satellite and 
linear multifilters. In the case of the kinematic SZ effect, 
to obtain V (v r ) it is necessary to divide the filtered map 
by y c , so error bars are larger for low comptonizations (as 
reflected by figure |^J. For clusters as considered here, with 
y c — 10~ 4 and temperature T e ~ 5 keV, the $ multifilters 
give statistical error bars ov — 0.26, that is, o~ Vr ~ 800 km 
s _1 . The situation becomes worse for fainter clusters. This 
means that Plank will not be able to tell us the velocities of 
individual clusters 1 . Since $ multifilters provide an unbiased 
estimator of the kinematic SZ effect it may be possible, how- 
ever, to determine mean peculiar velocities on large scales by 
averaging over many clusters. On the other hand, standard 
matched multifilters would lead to erroneous estimations of 
mean peculiar velocities due to their intrinsic bias. 

One of the strongest traits of linear filtering is that it re- 
quires only a small number of assumptions about the data to 
work. The multifilters need only three information elements 
to do their work: the frequency dependence of the SZ effects, 
a reasonable knowledge on the generic cluster profile (and 
the instrumental beam at each frequency) and the cross- 
power spectra of the different channels' background. Uncer- 
tainties and errors in the a priori information must be kept 
to a minimum and, when present, they should be taken into 
account. In the following we will discuss briefly some of the 
expected problematics that will arise in real life clusters, and 
suggest a few ways to deal with them. 



6.1 SZ frequency dependence 

The frequency dependence of the thermal and the kine- 
matic SZ effects are very well known, so few surprises will 
come from this direction. Second-order relativistic correc- 
tions, however, can play a role, changing the spectral shape 
of the effects. In particular, relativistic effects can change 
the crossover frequency of the thermal SZ effect, a region 
that is fundamental to the detection of the kinematic SZ 
effect. The sensitivity of Planck will be too low to constrain 
the optic al depth of relativistic electrons r re ; of individual 
clusters (lEnfilin fc Hansen! 1200411 . For other more sensitive 



At least if only multifilters are used to estimate the kinematic 
SZ effect. Maybe the combination of the multifiltcring technique - 
as a denoising step- with more sophisticated estimation methods 
already suggested in the literature could help to improve these 
results. This is an exciting open possibility for future work. 
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experiments it may be necessary to take this effect into ac- 
count. The relativistic correction to the thermal SZ effect 
frequency dependence is given bv |R ephaelil 1)19951) and it is 
a function of the frequency v and the temperature of the 
electrons in the intracluster gas T e . The relativistic correc- 
tion to the kinematic SZ effect is due to the Lorentz boost 
to the electrons caused by the bulk velocity, introduces very 
small spectral distortions on the kinematic SZ frequency 
dependence and for a 10 keV cluster moving at 1000 km 
s the effect is about an 8 % of the non-relativistic term. If 
T e is known, the relativistic corrections can be calculated, 
eqs. (1101 and U2t can be correspondingly modified and the 
method can be used without problems. 

If T e is not known, it can be considered as another 
quantity to be estimated. In principle, multifilters such as 
these her e described are able to estimate only amplitudes. 
However, Irlerranz et alj J2002bf) showed a way to estimate 
additional parameters by means of consecutive linear filter- 
ings. In that work, the additional parameter to be estimated 
was the core radius (as we will discuss later) and the filter 
under study was the scale- adaptive filter. It is straightfor- 
ward to generalise the idea to UMMF and other parameters 
such as the frequency dependence. The philosophy of the 
method is as follows: it can be shown that matched filters 
(an unbiased matched filters) give a maximum amplifica- 
tion of the sources when the correlation between the shape 
of the filter and the shape of the source is perfect. If for 
example the shape (profile) of the source is known except 
for a scale factor (such as the core radius), we can filter 
the image with a number of different matched filters with 
varying core radius parameter, and study how amplification 
changes with r c . The value of the test parameter at which 
the maximum amplification is obtained will be our estimate 
of r c . The same applies to multifilters: the maximum am- 
plification will be obtained when the frequency dependence 
assumed for the filter matches the true one. Then, we can 
pass a set of different multifilters with varying frequency de- 
pendences (in this case, parametrised by different values of 
T e ). The maximum amplification occurs when the value of 
T e of the multifilter coincides with the true one. 



that all the clu sters had the same cor e radius r c , which was 
known a priori. iHerranz et al.l (l2002bT) provided a method to 
deal with this uncertainty using the scale-adaptive filter on 
a single map. The generalizat ion to multifilters is straight- 
forward iHerranz et al.l2002af) . The basic idea is to filter the 
maps with a set of different filters varying the r c parameter 
and then to look for the maximum gain cluster by cluster. 
With this method it is possible to determine the value of r e 
with t ypical errors lower than the pixel size (IHerranz et alJ 
l2002bl) . 



6.2.2 Profile shape 

If instead of the modified multiquadric profile used in eq. 
1271 the cluster has a different profile the matching between 
the filter and the cluster will not be perfect and the filtering 
will loose performance. This will be an important issue in 
high resolution experiments where the cluster size is much 
larger than the beam size. In the case of Planck most clusters 
will be smaller than the beam size and after convolution 
their observed shapes will be dominated by the beam profile. 
Therefore, in the case of Planck the errors due to profile 
uncertainties will be small except for a few cases of very 
extense clusters. 

For these few extense cluster and for the case of high 
resolution experiments we suggest the use of an iterative 
technique to adaptively fit the true profile of a cluster in an 
empiric way. Starting with a given profile guess such as the 
one in eq. 1271 , the shape of the profile can be slowly varied 
in successive steps, introducing small filter shape perturba- 
tions and maximising the signal to noise gain after filtering. 
The maximum gain corresponds to the best fit between the 
profile of the filter and the profile of the cluster. A method 
for adaptive shape fitting that could be adapted to this prob- 
lem has been successfully accomplished in the modelling of 
gravitational len sing effect due to massi ve galaxy clusters 
(see for example iBroadhurst et al] J2004T) '). We will explore 
this possibility in a future work. 



6.2 Cluster profiles 

For the sake of clarity we have deliberately kept the exam- 
ples presented here as very simple, academic tests, at least 
regarding to our SZ cluster simulations. On the other hand, 
the background simulations -CMB, Galactic foregrounds, 
extragalactic point sources and Planck instrumental noise 
levels as well as pixel and beam sizes- are very realistic. 
The purpose of this is to focus on the effect of the "contami- 
nants" rather than the whole problem. Unfortunately, in real 
life clusters are not so nicely symmetric and we do not have 
a perfect a priori knowledge of their spatial profile. Luck- 
ily, we neither have a total ignorance about it. Normally we 
have some previous idea about how clusters look like. The 
main uncertainties here are: 



6.2.1 Size uncertainties 

Clusters of many different sizes do appear in CMB maps. 
Throughout the examples presented in this work we assumed 



6.2.3 N on- circular profiles 

In this work we have assumed symmetric circular profiles 
and beams. This condition can be relaxed without loss of 
generality. If the non-symmetric profile is previously known 
(for example thanks to high-resolution observations in X-ray 
or other wavelengths) specific non-symmetric filters can be 
calculated without problem. If the shape of the cluster is not 
previously known, an iterative scheme such as the previously 
mentioned can be attempted. 



6.2-4 Non iso-thermal clusters 

It is frequent to consider that T e is constant all along the 
cluster, which allows us to write nicely the kinematic SZ ef- 
fect as in eq. 11211 . But galaxy clusters are not iso-thermal. 
As a consequence of this, cluster parameters derived through 
the observation of SZ effect will be affected . iHansenl i2004bl) 
has shown that peculiar velocities will be systematically 
shifted by 10 — 20% if the inner structure of the cluster is 
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not considered in detail, which requires either observation- 
ally expensive X-ray observations or high angular resolution 
SZ observations. 



6.3 Background power spectrum 

The background power spectrum can usually be estimated 
directly from the data. This is a safe practise if SZ clusters' 
contribution to the total power spectrum is relatively small. 
Such is the case expected for the Planck mission. If the SZ 
clusters were so strong that they contributed significantly 
to the total power spectrum at all the scales, they would be 
so conspicuous that detection would not be a problem, and 
multifiltering would not be necessary. 



6.4 Extragalactic point source contamination 

Historically, extragalactic point sources have been the 
strongest source of contamination in SZ observations. In this 
work w e have included r ealist ic point source simulations us- 
ing the lToffolatti et alJ 11199a) model, whose predictions for 
the radio source counts have been rec ently confirmed by the 
WMAP mission JBennett et al.ll2003l) . We have shown that 
multifilter estimation is not substantially affected by these 
point sources. There are however two details that have not 
been included in the simulations and can affect the estima- 
tion of cluster parameters: the intrinsic galaxy clustering 
and the spatial correlation between individual galaxies and 
galaxy clusters. 



6.4-1 Point source clustering 

We used the point sources of the iToffolatti et alJ (J.1998) 
model assuming a uniform Poisson distribution of galaxies in 
the sky. In the real sky, however, galaxies are spatially corre- 
lated due to clustering. Temperature fluctuations due to un- 
resolved point sources are stronger if c lustering is cons i dered . 
Apart from increasing the error bars, lAghanim et al.l (I2004T) 
have shown that this effect leads to significant systematic er- 
rors in the determination of the thermal a nd kinematic S Z 
amplitudes by means of the SASZ method (Hansen 2 004af) . 
It is not clear how galaxy clustering will affect the multifil- 
tering method presented in this work. If the correlation scale 
of the galaxies is similar to the cluster size systematic shifts 
can occur 2 , but if the unresolved galaxies' fluctuation field 
is not correlated with the positions of the clusters the effect 
should be small. This is a problem worth to study in future 
works. 



6.4-2 Point sources associated to clusters 

Another factor not considered in the simulations is the as- 
sociations between clusters and individual galaxies. These 
associations can be due to the galaxies that form part of the 



2 This is due to the same kind of effect that is the topic of this 
work: a cluster-like structure (in this case a "lump" of corre- 
lated galaxies with a correlation scale approximately equal to the 
cluster size) is superimposed to the observed cluster, showing a 
different frequency dependence but nevertheless producing a bias 
in the estimation of the cluster parameters. 



cluster or to gravitational lensing, that increases the num- 
ber of faint galaxies that are observed in the direction of the 
cluster through magnification effects. Emission from such 
galaxies can fill in the SZE decrement, leading to a wrong 
estimation of the effect. Galaxies inside the cluster are likely 
to affect the estimation of SZ effect at high frequencies while 
gravitational lensing effect is expected to be an issue f or fre- 
quencies < 30 GHz llCarlstrom. Holder fc Reesell2002l) . The 
solution to this problem is increasing the resolution of the 
observations so that individual galaxies can be detected an 
accounted for; in low-resolution experiments such as Planck 
this will not be possible. The different frequency dependence 
of the SZ effect and galaxy emission will help to reduce con- 
tamination, but some residual effects will remain that can 
lead to additional unavoidable biases. 

We remark that the method presented in this work is 
not inconrpatibly with other SZ detection and estimation 
techniques. Linear filtering is an useful tool that can be used 
alone or as a step inside a more ambitious analysis of the 
data. The filters here introduced are computationally fast, 
robust, efficient and unbiased. They reduce the contamina- 
tion due to noise and other astrophysical emission and opti- 
mise the separation between the thermal and the kinematic 
effect using all the frequency channels available (and not 
only the channels around the thermal crossover frequency). 
In the context of the future Planck mission, the implemen- 
tation of the filters is straightforward, but due to the limi- 
tations on angular resolution and sensitivity the kinematic 
SZ effect will be badly determined. Experinrents with more 
sensitivity and angular resolution will be much better for 
the detection of this elusive effect, but the implementation 
of the filters will require much more care since the inner 
structure of the clusters will reveal its complexity. We have 
hinted in the discussion some ways to deal with this com- 
plexity; applications to specific cases will be addressed in 
future works. 
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